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ABSTRACT 

We report on eighteen months of multiwavelength observations of the blazar 3C 454.3 (Crazy Dia- 
mond) carried out in the period July 2007 - January 2009. In particular, we show the results of the 
AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008 - Jan- 
uary 2009. During the May 2008 - January 2009 period, the source average flux was highly variable, 
with a clear fading trend towards the end of the period, from an average 7-ray flux Fe>ioomoV ^ 
200 x 10 -8 photons cm~ 2 s _1 in May- June 2008, to i 7 E>iooMeV ~ 80 x 10 -8 photons cm -2 s -1 in 
October 2008 - January 2009. The average 7-ray spectrum between 100 MeV and 1 GeV can be fit 
by a simple power law, showing a moderate softening (from Tgrid ~ 2.0 to Tgrid ~ 2.2) towards the 
end of the observing campaign. Only 3 a upper limits can be derived in the 20-60 keV energy band 
with Super- AGILE, because the source was considerably off-axis during the whole time period. 

In July- August 2007 and May- June 2008, 3C 454.3 was monitored by RXTE. The RXTE/PCA light 
curve in the 3-20 keV energy band shows variability correlated with the 7-ray one. The RXTE/PCA 
average flux during the two time periods is i 7 3 - 20kcV = 8.4 x 10~ n erg cm~ 2 s _1 , and i 7 3_20keV = 
4.5 x 10 _11 erg cm~ 2 s _1 , respectively, while the spectrum (a power-law with photon index Tpca = 
1.65 ±0.02) does not show any significant variability. Consistent results are obtained with the analysis 
of the RXTE/HEXTE quasi-simultaneous data. 

We also carried out simultaneous Swift observations during all AGILE campaigns. Swift/XRT 
detected 3C 454.3 with an observed flux in the 2-10 keV energy band in the range (0.9 — 7.5) x 
10 -11 erg cm -2 s _1 and a photon index in the range Fxrt = 1-33 — 2.04. In the 15-150 keV energy 
band, when detected, the source has an average flux of about 5 mCrab. 

GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period in the radio, millimeter, 
near-IR, and optical bands. The observations show an extremely variable behavior at all frequencies, 
with flux peaks almost simultaneous with those at higher energies. A correlation analysis between the 
optical and the 7-ray fluxes shows that the 7-optical correlation occurs with a time lag of r — — 0.4+Qg 
days, consistent with previous findings for this source. 

An analysis of 15 GHz and 43 GHz VLBI core radio flux observations in the period 2007 July - 2009 
February shows an increasing trend of the core radio flux, anti-correlated with the higher frequency 
data, allowing us to derive the value of the source magnetic field. 

Finally, the modeling of the broad-band spectral energy distributions (SEDs) for the still unpub- 
lished data, and the behavior of the long-term light curves in different energy bands, allow us to 
compare the jet properties during different emission states, and to study the geometrical properties 
of the jet on a time-span longer than one year. 

Subject headings: galaxies: active - galaxies: quasars: general - galaxies: quasars: individual: 
3C 454.3 - galaxies: jets - radiation mechanism: non thermal 
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1. INTRODUCTION 

Among active galactic nuclei (AGNs), blazars show 
intense and variable 7-ray emission above 100 MeV 
(|Hartman et al.l I1999D . with variability timescales as 
short as a few days, or a few weeks. 

Blazars emit across several decades of energy, from 
the radio to the TeV energy band, and their spectral 
energy distributions (SEDs) typically show two distinct 
humps. The first peak occurs in the IR/Optical band in 
the Flat Spectrum Radio Quasars (FSRQs) and in the 
Low-energy peaked BL Lacs (LBLs), and at UV/X-rays 
in the High-energy peaked BL Lacs (HBLs). The sec- 
ond hump peaks at MeV-GeV and TeV energies in FS- 
RQs/LBLs and in HBLs, respectively. In the framework 
of leptonic models, the first peak is commonly interpreted 
as synchrotron radiation from high-energy electrons in 
a relativistic jet, while the second peak is interpreted 
as inverse Compton (IC) scattering of soft seed photons 
by the same relativistic electrons. A recent review of 
th e blazar emission mechanis ms and energetics is given 
in Celotti & Ghisellini (2008). Alternatively, the blazar 
SED can be explained in the framework of the hadronic 
models, where the relativistic protons in the jet are the 
primary accelerated particles, emitting 7-ray radiation 
by means of photo-pair and photo-pion pro duction (see 
iMucke fc Protheroeil200lt iMiicke et a"Dl2003l for a review 
on hadronic models). 

Since the launch of AGILE, the FSRQ 3C 454.3 
(PKS 2251+158; z = 0.859) became one of the most ac- 
tive sources in the 7-ray sky. Its very high 7-ray flux (well 
above 100 x 10 -8 photons cm" 2 s" 1 for E > 100 MeV), 
its flux variability (on a time scale of 1 or 2 days), and 
the fact that it was always detected during any AGILE 
pointing, made it earn the nickname of Crazy Diamond: 
3C 454.3 is no w playing the same role as 3C 279 had for 
EGRET (e.g., iHartman et"aHl2001al |bfl. 

Multiwavelength studies of variable 7-ray blazars are 
crucial in order to understand the physical processes 
responsible for the emission along the whole spec- 
trum. Sinc e the detection of the exceptional 2005 out- 



Gio mmi et all [20061 : iFuhrmann et~aTl |2006; 
2006), several monitoring campaigns were 



burst (see 
iPian et all 

carried out to follow the source multi frequency behay ior 
(|Viilata et al.l [2^061 [20071 IRaiteri et al.1 120071 l2008allbh . 
In mid July 2007, 3C 454.3 underwent a new opti- 
cal brightening, which triggered observations at all fre- 
quencies. AGILE performed a target of opportunity 
(ToO) re-pointing towards the source and detec ted it 
in a very high 7-ray state (jVercellone et al.ll2008L here- 
after V08). In November and December 2007, AGILE 
detected high 7-ray activity from 3C 454.3, triggering 
multiwav elength ToO campaigns , whose results are re- 
ported in IVercellone et all (|2009l . hereafter Paper I), in 



Donnarumma et al. | (|2009L hereafter Paper II), and in 
Anderhub et al.l (|2009f) . respectively. Paper I and Pa- 
per II demonstrated that to fit the simultaneous broad- 
band SEDs from radio to 7-ray data, inverse Comp- 
ton (IC) scattering of external photons from the broad 
line region (BLR) off the relativistic electrons in the 

jet was required . In an earlier work based on the 

IVercellone et alJ (l2007f) preliminary flu x estimate of the 
July 2007 flare, iGhisellini et al.l (|2007|) made a compari- 
son between the 3C 454.3 SEDs in 2000 (EGRET data), 



2005 (optical and X-ray flare), and 2007 (AGILE 7-ray 
flare), discussing the role of the bulk Lorentz factor T 
(associated with the emitting source compactness) dur- 
ing the different epochs. 

Moreover, the results of correlation analysis performed 
in Paper I was consistent with no time-lags between the 
7-ray and the optical flux variations. Such a result was 
recently confirmed by iBonning et al.l ([2009D who corre- 
lated optical, UV, X-ray and 7 -rav Q data. In a very 
recent paper, lAbdo et all (|2009l ) show the results of the 
first three months of Fermi observations of 3C 454.3, 
from 2008 July to October. They present for the first 
time the signature of a spectral break above a few GeV, 
interpreted as a possible break in the energy distribution 
of the emitting particles. 

In this paper (Paper III) we present both a re- 
analysis of the AGILE published data collected dur- 
ing the period July 2007 - December 2007, and the re- 
sults of multiwavelength campaigns on 3C 454.3 dur- 
ing a long-lasting 7-ray activity period between 2008 
May 10 and 2009 January 12. In particular, we show 
the results of the AGILE campaigns which took place 
on May- June 2008 (mj08), July- August 2008 (ja08), 
and October 2008 - January 2009 (oj09). Preliminary 
7-ray r esults were distribut e d in Donnarumma et ahl 
([200811: iVittorini et al.1 (poM ): iGasparrini et al.l (|2008j k 
iPittori et al.l (l2008f). whi l e rad io-to-optical data were 
published in lVillata et"aH (|2009[ ). 

The paper is organized as follows: in Sect. [2] through [6] 
we present the AGILE, Swift, Rossi X-ray Timing 
Explorer (RXTE), GL AST- AGILE Support Program 
within the Whole Earth Blazar Telescope (GASP- 
WEBT) and radio VLBI data analysis and results; in 
Sect. [7] we present the simultaneous multiwavelength 
light-curves. In Sect. [5] and M we discuss the results and 
draw our conclusions. Throughout this paper the quoted 
uncertainties are given at the la level, unless otherwise 
stated, and we adopted a A-CDM cosmology with the fol- 
lowing values for the cosmological parameters: h = 0.71, 
fl m = 0.27, and fl A = 0.73. 

2. AGILE DATA 
2.1. Data Reduction and Analysis 

The AGILE satellite (|Tavani et al.ll2008l [20091) . a mis- 
sion of the Italian Space Agency (ASI) devoted to high- 
energy astrophysics, is currently the only space mission 
capable of observing cosmic sources simultaneously in the 
energy bands 18-60 keV and 30 MeV-50 GeV. 

The AGILE scientific instrument combines four ac- 
tive detectors yielding broad-band coverage f rom hard X- 
rays to gamma-rays: a Silicon Tracker (ST; iPrest et ahl 
12001 30 MeV-50 GeV), a co-ali gned coded-mask hard 
X-ray imager, Super- AGILE (SA: lFeroci eTaD[2007l 18- 
60 keV), a non- ima ging Csl Mini-Calorimeter (MCAL; 
lLabanti et all [200l 0.3-100 Me V), and a segmented 
Anti- Coincidence System (ACS; iPerotti et all 120061) . 
Gamma-ray detection is obtained by the combination of 
ST, MCAL and ACS; these three detectors form the AG- 
ILE Gamma- Ray Imaging Detector (GRID). 

53 The 7-ray data are taken from the Fermi/LAT 
monitored source list light curves available at 
http : //f ermi . gsfc .nasa. gov/ssc/dat a/access/ 
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TABLE 1 
AGILE/GRID observation log. 



Level- 1 AGILE-GRID data were analyzed using the 
AGILE Standard Analysis Pipeline (see V08 for a de- 
tailed description of the AGILE data reduction). Since 
3C 454.3 was observed at high off-axis angle due to the 
satellite solar panel constraints, an ad-hoc 7-ray analysis 
was performed. We used 7-ray events filtered by means 
of the FM3.119_2a AGILE Filter pipeline. Counts, ex- 
posure, and Galactic background 7-ray maps were cre- 
ated with a bin-size of 0?25 x 0?25, for E > 100 MeV. 
Since the source was observed up to 40° off-axis, all the 
maps were generated including all events collected up to 
60° off-axis. We rejected all 7-ray events whose recon- 
structed directions form angles with the satellite-Earth 
vector smaller than 85°, reducing the 7-ray Earth albedo 
contamination by excluding regions within ~ 15° from 
the Earth limb. We used the version (BUILD- 16) of the 
Calibration files (10006), which are publicly available at 
the ASI Science Data Centre (ASDC) siteg) and of the 
7-ray diffuse emission model (| Giuliani et al. I [20041 ). We 
ran the AGILE Source Location task in order to derive 
the most accurate location of the source. Then, we ran 
the AGILE Maximum Likelihood Analysis (ALIKE) using 
a radius of analysis of 10°, and the best guess position 
derived in the first step. The particular choice of the 
analysis radius is dictated to avoid any possible contam- 
ination from very off-axis residual particle events. 



2.2. GRID Results 

Table Q] shows the AGILE/GRID observation log dur- 
ing the different time periods. We have re-analysed all 
the AGILE data already published in V08, Paper I, and 
Paper-II, respectively, according to the procedure de- 
scribed in Section 12.11 The results are discussed in Sec- 
tion [7j In the following paragraphs we report the 
detailed results of the still unpublished AGILE data. 

2.2.1. May- June 2008 

The AGILE campaign was split into two different pe- 
riods, May 10-June 9 (PI) and June 15-30 (P2) be- 
cause of a ToO re-pointing towards W Comae. The 
total on-source exposure is 1.03 (PI) + 0.54 (P2) Ms. 
3C 454.3 was detected, during PI and P2, at a 25.6ct 
and 16.3<7 level with an average flux of -Pe>ioomcV = 
(218 ± 12) x 10" 8 photons cm" 2 s~\ and ^|>i 00 McV = 
(198 ±17) x 10 _8 photons cm~ 2 s _1 , respectively, as de- 
rived from the AGILE Maximum Likelihood Code anal- 
ysis. 

54 http://agile.asdc.asi.it 



Figure [TJ filled circles in panel (a), shows the 7- 
ray light-curve at 1-day resolution for photons above 
100 MeV. 
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Fig. 1— Panel (a): AGILE/GRID 7-ray light curve at « 1- 
day resolution for E>100 MeV in units of 10 — 8 photons cm -2 s — 1 
during the period May 2008 - January 2009. The downward arrows 
represent 2cr upper-limits. Panel (b): the same as in panel (a), but 
with a time bin of three (filled circles and squares), and seven (filled 
triangles) days, respectively. 



We note that, particularly at the beginning of the cam- 
paign, 3C 454.3 was not always detected on a day-by-day 
timescale. On MJD ~ 54610 the source began to be 
detected at a 3cr level almost continuously; this clearly 
indicates the onset of a 7-ray flaring activity. 

The average 7-ray flux as well as the daily values were 
derived a ccording to the 7-ra y analysis procedure de- 
scribed in iMattox et al.l (|1993l ). First, the entire period 
was analyzed to determine the diffuse gas parameters and 
then the source flux density was estimated independently 
for each of the eighteen 1-day periods with the diffuse pa- 
rameters fixed at the values previously obtained. 

Figure Q] panel (b), shows the same AGILE/GRID 
data binned on a time scale of three days. The light 
curve clearly shows a strong degree of variability, with a 
dynamic range of about four in about two weeks. 

Figure [21 panel (a) , shows the average 7-ray spectra 
extracted over the observing periods PI and P2. 

Each average spectrum was obtained by computing the 
7-ray flux in five energy bins over each period: 50 < E < 
100 MeV, 100 < E < 200 MeV, 200 < E < 400 MeV, 
400 < E < 1000 MeV, and 1000 < E < 3000 MeV. We 
note that the current instrument response is accurately 
calibrated in the energy band 100 MeV-1 GeV, and that 
the flux above 1 GeV is underestimated by a factor of 
about 2. For those reasons, we fit the data by means 
of a simple power law model and restricted our fit to 
the 100 MeV-1 GeV energy range, obtaining (in units of 
photons cm~ 2 s _1 MeV -1 ): 



F P1 (E) = 2.63 x 10" 



F^(E) = 1.58 x 10" 



E 



I MeV") 



(2.05±0.10) 



E 

lMeV 



-(1.98±0.16) 



(1) 



(2) 
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Fig. 2. — Panel (a): AGILE/GRID average 7-ray spectrum for 
periods PI and P2. The blue-dashed and the red-dotted lines 
represent the best-fit power law models for PI and P2, respec- 
tively. Panel (b) and Panel (c) show the average 7-ray spectra 
during the periods July-August 2008 and October 2008 - Jan- 
uary 2009, respectively. In the three panels only three energy bins 
were considered for the spectral fitting: 100 < E < 200 MeV, 
200 < E < 400 MeV, 400 < E < 1000 MeV (see text for details). 
[See the electronic edition of the Journal for a color version of this 
figure.) 



The different energy range, and, above all, the differ- 
ent time period, could explain the different value of the 
AGILE (2.05 ± 0.10, a nd 1.98 ± 0.16) an d Fermi-LAT 
(2.27 ± 0.03 , pre-break, lAbdo et al.ll2009D 7-ray photon 
indices. 

2.2.2. July- August 2008 

The AGILE campaign started immediately after the 
Fermi/LAT detection of a ve ry high 7-ray act ivity in 
the period 2008 July 10-21, (fTosti et alj|2008l ). which 
reached, on July 10, a 7-ray flux of ^E> r iu0MoV = 
1200 x 10" 8 photons cm" 2 s" 1 (jAbdo et all 12009}). The 
AGILE observations covered the period from 2008-07- 
25 19:57 UT to 2008-08-14 21:08 UT, for a total on- 
source exposure of about 0.71 Ms. 3C 454.3 was detected 
at a 17.5c level with an average flux of -Fe>ioomcV = 
(255 ±21) x 10 -8 photons cm -2 s _1 , as derived from the 
AGILE Maximum Likelihood Code analysis. 

Figure [l] filled squares in panel (a) and in panel (b) , 
shows the 7-ray light-curve at 1-day and at 3-day res- 
olution, respectively, for photons above 100 MeV. The 
average 7-ray flux as well as the daily values were de- 
rived according to the procedure described in § 12.2.11 . 
Contrary to the May- June data, the July-August light 
curve does not show any clear sign of variability. 

Figure [21 panel (b), shows the average 7-ray spectrum 
derived over the entire observing period. The average 
spectrum was obtained by computing the 7-ray flux in 
the same way as in § 12.2.11 We fit the data by means of 
a simple power law model and restricted our fit to the 



most reliable energy range (100 MeV-1 GeV): 

g \ -(2.11±0.14) 



F }MS {E) = 3.96 x 10" 



MeV/ 



, (3) 



in units of photons cm~ 2 s _1 MeV -1 . During this pe- 
riod, which partially overlaps with the Fermi one, the 
AGILE 7-ray photon index is, within the statistical er- 
rors, in agreement with the Fermi-LAT result. 

2.2.3. October 2008 - January 2009 

The AGILE observations covered the period from 2008- 
10-17 12:51 UT to 2009-01-12 14:30 UT, for a total on- 
source exposure of about 2.86 Ms. 3C 454.3 was detected 
at a 17.9cr level with an average flux of -Fe J ^'j' 00McV = 

(77 ± 5) x 10 -8 photons cm -2 s -1 , as derived from the 
AGILE Maximum Likelihood Code analysis. 

Figure [TJ filled triangles in panel (a), shows the 7- 
ray light-curve at 1-day resolution for photons above 
100 MeV. The average 7-ray flux as well as the daily val- 
ues were derived according to the procedure described in 
§ 12.2.11 The light curve does not show any clear trend, 
partly because of a dominant fraction of upper limits in 
the data. For this reason, we decided to rebin the light 
curve on a time scale of one week. The resulting light 
curve is shown in Figure [TJ filled triangles in panel (b) . 
On this time scale, a clear trend is present, with the 
source dimming as a function of time, with a dynamic 
range of about a factor of two. 

Figure [H panel (c), shows the average 7-ray spectrum 
derived over the entire observing period. The average 
spectrum was obtained by computing the 7-ray flux in 
the same way as in § 12.2.11 We fit the data by means of 
a simple power law model and restricted our fit to the 
most reliable energy range (100 MeV-1 GeV): 



F° m (E) = 2.10 x 10" 



t— , \ -(2.21±0.13) 

) , (4) 



lMeV 



in units of photons cm 2 s 1 MeV 1 

2.3. Super- A GILE Results 

During the various AGILE pointings, 3C 454.3 was lo- 
cated substantially off-axis in the Super-AGILE field of 
view (FoV). For this reason, only 3c upper limits can 
be derived in the 20-60 keV energy band during the AG- 
ILE/GRID observations. Table [2] summarizes the Super- 
AGILE observations and their results. 

3. Swift DATA 

3.1. Data Reduction and Analysis 

Swift pointed observations (jGehrels et al.l I2004D were 
performed from 2007-07-26 to 2009-01-01. These obser- 
vations were obtained both by means of several dedicated 
ToOs (PI S. Vercellone) and by activating Swift Cycle-3 
(Obs. ID 00031018 PI A.W. Chen) and Cycle-4 Propos- 
als (Obs. ID 00031216, PI S. Vercellone). A long-lasting 
monitoring program (P.Is. L.Fuhrmann and S. Vercel- 
lone) covers the period July-October 2008. 

3.1.1. Swift/ XRT 

Table [3] summarizes the Swift/XRT observations. The 
XRT data were processed with standard procedures 
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TABLE 2 

SUPER-AGILE OBSERVATION RESULTS. 



Start Time End Time #x $z Exposure i"20-60keV 

(UTC) (UTC) (Deg.) (Deg.) (ks) (mCrab) 

2008-05-31 10:18 2008-06-09 13:38 -23.0 +06.0 380 < 16 

2008-06-15 14:11 2008-06-21 12:59 -36.0 +08.0 270 < 18 

2008-07-25 21:39 2008-08-02 23:29 +03.4 -42.0 345 < 18 

2008-10-17 18:47 2008-10-29 23:12 -00.8 -45.0 460 < 21 



(xrtpipeline vO.12.1), adopting the standard filter- 
ing and screening criteria, and using FTOOLS in the 
Heasoft package (v. 6. 6.1). 

The source count rate was variable during the cam- 
paigns, ranging from 0.26 to 1.8 counts s _1 . For this 
reason, we considered both photon counting (PC) and 
windowed timing (WT) data, and further selected XRT 
event grad es 0-12 and 0-2 for th e PC and WT events, re- 
spectively (jBurrows et al.ll2005fl . Several Swift/XRT ob- 
servations showed an average count rate of > 0.5 counts 
s _1 , therefore in these cases pile- up correction was re- 
quired for the PC data. We extracted the source events 
from an annular region with an inner radius of 3 pixels 
(estimated by means of the PSF fitting technique) and 
an outer radius of 30 pixels (1 pixel ~ 2'.'36). When the 
average count rate was < 0.5 counts s _1 , we used the full 
30-pixel radius region. 

We also extracted background events within an annu- 
lar region centered on the source with radii of 110 and 
116 pixels. Ancillary response files were generated with 
xrtmkarf , and account for different extraction regions, 
vignetting and PSF corrections. We used the spectral re- 
distribution matrices vOll in the Calibration Database 
maintained by HEASARC. Swift/XRT uncertainties are 
given at 90% confidence level for one interesting param- 
eter (i.e., Ax 2 = 2.71) unless otherwise stated. 

The Swift/XRT spectra were rebinned in order to have 
at least 20 counts per energy bin. We fit the spectra 
with an absorbed power law model, (wabs* (powerlaw) 
in XSPEC 11.3.2). The Galactic absorption was fixed to 
the value of i Vff al = 1.34 x 10 21 cm~ 2 , as obtained by 
IVillata et all (l200ll by means of a deep Chandra obser- 
vation. We note the adopted value is consistent with the 
mean of the distribution of the Nn values obtained by 
fitting the spectra with an absorbed power law model 
and free absorption. 

3.1.2. Swift/UVOT 

The UVOT data analysis was performed using the 
uvotimsum and uvotsource tasks included in the FTOOLS 
software package (HEASOFT v6.6.1). The latter task 
calculates the magnitudes through aperture photometry 
within a circular region and applies specific corrections 
due to the detector characteristics. Source counts were 
extracted from a circular region with a 5 arcsec radius. 
The background was extracted from source-free circular 
regions in the source surroundings. The reported mag- 
ni tudes are on the U VOT photometric system described 
in iPoole et"aT1 (|2008| ) , and are not corrected for Galactic 
extinction. 



We analyzed Swift/RAT Survey data in order to study 
the hard X-ray emission of 3C 454.3 and to investigate 
its evolution as a function of time. 

We produced a light curve for the sou rce at a 16-d 
binning usin g the procedures described in Krim m~et al.l 
(|2006t I2008L and references therein; also seeE3). 



3.2. Results 

Figure [3] shows the Swift/XRT fluxes in the 2-10 keV 
energy range, the Swift/XJVOT observed magnitudes 
(in the V, B, U, Wl, M2, and W2 bands), and the 
Swift/RAT fluxes in the 15-150 keV energy range as a 
function of time for the whole observing period. 
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3.1.3. Swift/ 'BAT 



http : //swift . gsf c .nasa. gov/docs/swif t /result s/transients/Transieni 
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TABLE 3 

Swift/XKT OBSERVATION LOG. 



Sequence 


Obs. 


Start time (UT) 


End time (UT) 


Exp. a 




Mode 


(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


s 


(1) 


(2) 


(3) 


(4) 


(5) 


00035030013 


PC 


2007-07-26 00:55:44 


2007-07-26 01:13:58 


1073 


00035030014 


PC 


2007-07-28 07:26:46 


2007-07-28 10:44:55 


817 


00035030015 


PC 


2007-07-30 10:59:09 


2007-07-30 14:16:56 


897 


00035030016 


PC 


2007-08-01 11:05:10 


2007-08-01 11:07:58 


168 


00035030017 


WT 


2007-08-01 09:33:17 


2007-08-01 13:06:59 


3903 



1 NOTE.-Table[3]is published in its entirety in the electronic edition of 
the Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 

a The exposure time is spread over several snapshots during each obser- 
vations 



Fig. 3. — Panel (a): Swift /UVOT light curves (observed magni- 
tudes) in the V (red triangles), B (green quares), and U (blue cir- 
cles). Panel (b): Swift/WOT light curves (observed magnitudes) 
in the Wl (red triangles), M2 (green quares), and W2 (blue cir- 
cles). Panel (c): Swift/XRT light curve (observed fluxes) in the 
2-10 keV energy band and in units of 10 — 11 erg cm" 2 s -1 . Panel 
(d): Swift/BAT light curve in units of mCrab in the energy band 
15 — 150 keV. Downward arrows mark 3<r upper limits. 



In order to diminish the statistical uncertainties, we se- 
lected observations with a number of degrees of freedom 
(dof) > 10. We note that a common dimming trend is 
present both in the UV and in the X-ray energy bands. 

As shown in Figure [31 panel (d) , the source has not 
been always detectable by Swift/BAT throughout the 
considered period, and in several time interval only 3<r 
upper limits can be derived. 

Figure |4] shows the Swift/XRT photon index as a func- 
tion of the X-ray flux in the 2-10 keV energy band. Black 
circles and red squares represent data acquired in PC and 
WT mode, respectively. WT data are not affected by 
pile- up at the observed count rate (CR < 3 counts s -1 ). 
We investigated the possible presence of a spectral trend 
in the X-ray data. If we consider WT data only, a 
"harder- when-brighter" trend seems to be present. Fit- 
ting the data with a constant model, we can exclude this 
model at the 99.9993% level. When analyzing the PC 
data only (as well as the sum of the PC and WT data), 
this spectral trend vanishes, and a fit with a constant 
model still holds. Nevertheless, if we exclude the points 
at fluxes F2_io kov < 2 x 10 -11 erg cm -2 s _1 a trend still 
holds. These points at low fluxes could correspond to 
physically different state of the source than the high 
fluxes one. 

We also note that a deep and prolonged monitoring 
of 3C 454.3 at mid and low X-ray states (-Fi-iokeV J$ 
10~ n erg cm~ 2 s" 1 ) will be crucial to test the possible 
presence of a spectral trend. Our data-set contains only 
four observations (90023002, 90023003, 90023006, and 
90023008) at this flux level, which were acquired during 
the source low state in December 2008. 

4. RXTE DATA 

4.1. Data Reduction and Analysis 

The Rossi X-ray Timing Explorer (RXTE) satellite 
observed 3C 454.3 in two epochs: from 2007-07-28 to 
2007-08-04 and from 2008-05-30 to 2008-06-19. Here we 
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Fig. 4. — Swift/XRT photon index as a function of the 2-10 keV 
flux. Red squares and black circles mark the Swift/XRT windowed 
timing (WT) and photon counting (PC) data, respectively. 



report the analysis of the data obt ained both with th e 
Proportional Counter Array (PCA, IJahoda et al.|[l996D . 
which is sensitive in the 2-60 keV energy range, and with 
the High-Energy X-Ra y Timing Experiment (HEXTE, 
iRothschild et al.l Il998i ). which is sensitive in the 15- 
250 keV energy range. RXTE data were collected by 
activating a Cycle 12 ToO proposal (ID. 93150, PI 
AW. Chen). 

The PCA is composed of 5 identical units (Proportional 
Counter Units, PCUs), but during our observations only 
part of them were used. Since PCU2 was the only unit 
always on during our observations and it is the one which 
is best calibrated, we report the results obtained from 
the PCU2 data only. The data were processed using 
the FTOOLS v6.4.1 and screened using standard filtering 
criteria. The net exposure times for the whole data-set 
in the first and second epoch were 36.6 ks and 17.4 ks, 
respectively. 

The background light curves and spectra for each ob- 
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servation were produced using the model appropriate to 
faint sources. We restricted our analysis to the 3-20 keV 
energy range, in order to minimize the systematic er- 
rors due to background subtraction and calibration of 
the PCA instrument. 

Figure [5] shows the 3-20 keV light curve of the whole 
RXTE/PCA data-set. 
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Fig. 5.— Panel (a) and (b) show the RXTE/PCA light curve 
in the energy band 3-20 keV during the periods 2007 July 29 - 
August 5 (black points) and 2008 May 30 - June 19 (red triangles), 
respectively. [See the electronic edition of the Journal for a color 
version of this figure.] 

Strong variability is observed when comparing the 
count rates of different observations. Moreover, the av- 
erage count rate during the second epoch (panel (b)) is 
reduced. In order to investigate possible changes in the 
spectral shape with time we extracted light curves in two 
energy ranges (3-7 keV and 7-20 keV). Their hardness 
ratio did not show any significant variation. 

A cumulative spectrum for the first and the second 
epoch was extracted and simultaneously fitted with a 
power-law model corrected for photoelectric absorption 
(wabs* (powerlaw) model in XSPEC), allowing only the 
power-law normalization to assume a different value in 
the two spectra. Figure [6] shows the RXTE/PCA spec- 
tra for both periods. A good fit (x 2 = 68.3 for 76 degrees 
of freedom) was obtained with the following best-fit pa- 
rameters (errors are at the 90% confidence level): pho- 
ton index T = 1.65±0.02, and a flux in the 3-20 keV 
energy band i 7 3„ 2 okcV = 8.4 x 10~ n erg cm~ 2 s _1 and 
-F3-20 keV = 4.5 x 10~ n erg cm~ 2 s _1 for the first and 
second epoch spectrum, respectively. 

We note that the average RXTE/PCA flux during the 
7-ray flare detected in July 2007 was about a factor of 
two higher than the flux detected about 10 months later. 
Moreover, during both the July 2007 and the May- June 
2008 campaigns, the hard X-ray flux varied significantly, 
by about 50%, on a time scale of about one week. 

We also analyzed the data of the HEXTE. Only the 
data from cluster B were analyzed, since the rocking sys- 
tem for background evaluation was disabled in the other 
instrument cluster. After a standard processing^ we ex- 
tracted an average spectrum from all the available data, 

56 http: //heasarc .gsf c .nasa. gov/docs/xte/recipes/hexte .html 
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Fig. 6.— RXTE/PCA average spectra for both periods, July 2007 
(black points), and May- June 2008 (red triangles), respectively. 
[See the electronic edition of the Journal for a color version of 
this figure.] 

for a dead-time corrected exposure time of 18 ks. The 
source was detected up to ~ 50 keV and its spectrum can 
be fit well (x 2 =15.1 for 19 d.o.f.) by a power-law model 
with a photon index of 1.6 ± 0.1, in perfect agreement 
with the photon index derived from the PCA spectrum. 
Also, the normalization of the HEXTE spectrum is con- 
sistent with an high energy extrapolation of the time- 
averaged PCA spectrum: the observed flux in the 20- 
40 keV energy band is (5±3)xl0 -11 erg cm~ 2 s _1 . This 
flux (approximately 6 mCrab) is also consistent with the 
upper limits obtained by Super-AGILE in the same time 
periods. 

5. OPTICAL-TO-RADIO DATA 
5.1. GASP-WEBT Data Reduction and Analysis 

The Whole Earth Blazar Telescope (WEBT0 has 
been monitorin g 3C 454.3 since the exceptional 2004- 
2005 outbur st (|Villata et all 120061 l2007t iRaiteri et all 
120071 l2008a lTbh. throughout the whole period of the A G- 
ILE ob servation. We refe r to IRaiteri et all (|2008allbl ) 
and to iVillata et al.l ((2009) for a detailed presentation 
and discussion of the radio, mm, near-IR, optical and 
Swift/UVOT data. 

Figure [7J shows the GASP-WEBT light curve in the 
R optical band, displaying several intense flares with a 
dynamic range of ~ 2.4 mag in about 14 days, while Fig- 
ure[8]shows the GASP-WEBT light curves in the near-IR 
(J, H, K), radio (5, 8, and 14.5 GHz), and mm (37, 230, 
and 345 GHz), respectively^! 

6. RADIO VLBI DATA 

6.1. Radio VLBI Data Reduction and Analysis 

High resolution radio VLBI data were obtained from 
the MO JAVE (Monitoring Of Jets in Active galactic nu- 
clei with VLB A Experiments) project, a long-term pro- 
gram to monitor radio brightness and polarization varia- 
tions in jets associated with active galaxies visible in the 

57 http : //www , oato ■ inaf . it /blazar /webt, see e.g. 
IVillata et atlpOQl ). 

58 The radio-to-optical data presented in this paper are stored 
in the GASP-WEBT archive; for questions regarding their avail- 
ability, please contact the WEBT President Massimo Villata 
(villataSoato . inaf . it). 
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Fig. 7.— GASP-WEBT light curve in the R optical band in the 
2007-2008 and 2008-2009 observing seasons. 



16 p 1 1 1 , 1 1 , 1 , 1 n 




54400 54600 
Time [MJD] 



Fig. 8. — Panel (a): low-frequency radio data. Red triangles, 
blue squares and black circles represent the radio flux at 5, 8, and 
14.5 GHz, respectively. Panel (b): high-frequency radio data. Red 
triangles, blue squares, and black circles represent the radio flux at 
37, 230, and 345 GHz, respectively. Panel (c): near-IR data. Red 
triangles, blue squares, and black circles, represent the J, H, and 
K bands, respectively [See the electronic edition of the Journal for 
a color version of this figure.] 

northern sky (|Lister et al.1 120091 : see alscPl). The object 
was observed with the full VLB A at 15 GHz. We ob- 

59 http : //www . physics . purdue . edu/MO J AVE 



tained the calibrated / images and used the MPS pack- 
age to derive the position and flux density of the core and 
of a few substructures in the jets using the task JMFIT 
(Gaussian fit) (see Figure [9|). 
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Fig. 9.— VLBI image of 3C 454.3 at 15 GHz on 2007 August 9 
(MJD 54321). The peak flux density is 2.8 Jy beam -1 and con- 
tours are traced at ±(1, \/{2), 2, ...) X 1.0 mjy beam — 1 . The 
cross in the bottom left corner shows the beam FWHM, which 
is 1.07 X 0.52 mas at —5.4 deg. 

Moreover, this source was additionally observed by 
VLBA at four epochs during the period of the maxi- 
mum brightness within the BK150 VLBA experiment 
to measure parsec-scale spectra of 7-ray bright blazars 
(Sokolovsky et al., in preparation). We use 15 and 
43 GHz results from this program to provide better ra- 
dio coverage of the high activity period. These data are 
in agreement with MOJAVE results and give a better 
statistics in the high active period. 

The core is always unresolved by our Gaussian fit and 
uncertainties on the flux density are dominated by cali- 
bration uncertainties (a few percent). 

In Figure |TD] we show the 3C 454.3 VLBI radio core 
flux (panel (a)) at 15 and 43 GHz, the radio components 
flux density at 15 GHz (panel (b)), and the distance of 
radio components from the core (panel (c)) as a function 
of time. 

The flux shows a constant increase from 2006 June 15 
(MJD 53901) till 2008 October 3 (MJD 54742), followed 
by a fast decrease towards the last epoch presented here, 
2009 June 25 (MJD 55007). Jet components show a well 
defined flux density decrease (component 1) or a slower 
flux density decrease which becomes almost constant in 
the last epochs. Proper motion is evident, but slowing 
in time for components 1 and 2; it is almost absent for 
component 3. 

All data are in agreement with a strong core flux den- 
sity variability possibly connected to the 7-ray activity, 
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( 2009), and iDonnarumma et al.l ()2009l ) . respectively. 




3C 454.3 Jul. 2007 - Jan. 2009 
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Fig. 10. — Panel (a): radio core flux density at 15 GHz (filled cir- 
cles) and at 43 GHz (open squares), respectively. Panel(b): radio 
components flux density at 15 GHz. Panel (c): radio components 
motion at 15 GHz. The vertical dashed lines represent the start 
(2007 July 24) and the stop (2009 January 12) of all AGILE ob- 
servations, respectively. 



while jet components are moving away and slowly de- 
creasing in flux density, and are not affected by the re cent 
core activity. In recent paper, IKovalev et al.l (|2009D in- 
deed find a connection between the radio and the 7-ray 
emission, correlating the Fermi three month data with 
the MO J AVE ones, and arguing that the central region 
of the blazars being the source of 7-ray flares. Neverthe- 
less, a detailed study of the radio structure of 3C 454.3 
is beyond the aims of this paper, therefore the jet prop- 
erties will be discussed in depth elsewhere (Lister et al. 
2009 in preparation, and Jorstad et al. 2010 in prepara- 
tion) . For this reason in the following we will concentrate 
only on the core. 

In the last two years this source has also been observed 
with the VLBA at 43 GHz (Jorstad et al. 2010; see also 
FT ) . We used the available images to derive the core flux 
density of the core at 43 GHz. Note that, for a better 
comparison with 15 GHz VLBI data, at 43 GHz we used 
natural weights and we have not searched for possible 
core subcomponents (we refer to Jorstad et al. 2010 for 
a detailed study of the radio structure) . 

The radio core shows an inverted spectrum (self- 
absorbed), more evident in the high active regime, fol- 
lowed by a strong flux density decrease. In this region 
the radio spectrum is no longer inverted. 

7. EIGHTEEN MONTHS OF MONITORING 

In this section we present a summary of all the ob- 
servations on 3C 454.3 in the period between 2007 July 
24 and 2009 January 12. The results of the campaigns 
performed in 2007 July, Nove m ber a n d December were 
discussed in iVercellone et al.l (|2008fl . iVercellone et al.l 

60 http : //www. bu. edu/blazars/VLBAproject .html 
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Fig. 11.— AGILE/GRID light curve at fa 3-day resolution for 
E>100 MeV in units of 10 -8 photons cm s . Different colors 
correspond to different observing campaigns, as described in Ta- 
ble [0 

Figure [TT] shows the AGILE/GRID light curve at 
~ 3-day resolution for E>100 MeV in units of 
10 -8 photons cm~ 2 s _1 . The light curve shows several 
7-ray flares, with a dynamical range of a factor of 3-4 
on a time scale of about ten days. Moreover, a clear 
dimming trend in the long-term light curve is present. 
Table [4] shows the AGILE/GRID fluxes and spectral in- 
dices derived at different epochs. 

7.1. Multiwavelength light curves 

The AGILE/GRID wide field of view allowed for the 
first time a long-term monitoring of 3C 454.3 at ener- 
gies above 100 MeV. Moreover, coordinated and almost 
simultaneous GASP-WEBT and Swift observations pro- 
vided invaluable information on the flux and spectral be- 
havior from radio to X-rays. 

Figure [T2l shows the 3C 454.3 light curves at different 
energies over the whole period. 

The different panels show, from bottom to top, the 
AGILE/GRID light curve at ~ 1-day resolution for 
E>100 MeV in units of 10~ 8 photons cm~ 2 s -1 , the 
Swift/HAT light curve in the energy range 15-150 keV at 
« 2-weck resolution, the Swift/XBT (filled circles) and 
the RXTE/PCA (filled squares) light curves in the en- 
ergy range 3-10 keV , the Swift/WOT light curve in the 
UV W2 filter, the Swift/WOT light curve in the optical 
B filter, the GASP-WEBT light curve in the optical R 
filter, and the VLBI radio core at 15 GHz (filled circles) 
and the UMRAO 14.5 GHz (open circles) light curves, 
respectively. 

We note that RXTE/PCA data are systematically 
higher than Swift/XRT ones, which is consistent with the 
20% uncertainty in the relative calibrations of the two in- 
strum ents in this energy band, reported bv lKirsch et al.l 
(j20p. 

Figure [T2] shows the light curves in the R band, at 
1.3mm (230 GHz), and above 100 MeV. The light curve 
in the millimeter wavelength shows a different behav- 
ior starting from the enhanced 7-ray activity at MJD^ 
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TABLE 4 

AGILE/GRID 7-ray fluxes and spectral indices above 100 MeV at 

DIFFERENT EPOCHS. 



Start Time 


End Time 


Fe>100 McV 


r 


(UTC) 


(UTC) 


(xlO -8 phcrn" 2 s -1 ) 





2007-07-24 14:30 2007-07-30 11:40 416.2 ± 36.0 1.74 ±0.16 

2007-11-10 12:16 2007-12-01 11:38 224.2 ± 15.3 1.91 ±0.14 

2007- 12-01 11:39 2007-12-16 12:09 265.7 ±17.5 1.86 ± 0.12 

2008- 05-10 11:00 2008-06-09 15:20 218.5 ± 12.2 2.05 ± 0.10 
2008-06-15 10:46 2008-06-30 11:14 198.5 ± 17.1 1.98 ±0.16 
2008-07-25 19:57 2008-08-14 21:08 254.8 ± 20.6 2.11 ± 0.14 
2008-10-17 12:51 2009-01-12 14:30 77.0 ± 5.5 2.21 ± 0.13 



54600, as will be discussed in Sect. 18.51 

Moreover, starting from MJD 54750, the whole 
jet seems to become less energetic, with an almost 
monotonic flux decrease, except for a minor burst at 
MJD 54800. 

8. DISCUSSION 

In the following Sections, we will discuss the corre- 
lations between the flux variations in different energy 
bands, the properties of the jet, and the physical pa- 
rameters of the emitting source. This latter point will 
be addressed by means of complementary approaches, 
namely the SED model fitting and discussing the geo- 
metrical properties of the jet itself. 

8.1. Correlation analysis 

We investigated the correlation between the 7-ray flux 
and the optical flux density in the R band by means of 
the d i screte correlation function ( DCF. lEdelson fc Krolikl 
Il988t iHufnagel fc Bregmanlll992t ). This method was de- 
veloped to study unevenly sampled data sets and can 
give an estimate of the accuracy of its results. Because 
of the sampling gaps in the light curves, especially at 
7-rays, we calculated the DCF on four distinct periods: 
July 2007 (mid 2007), November-December 2007 (Fall 
2007), May-August 2008 (mid 2008), and October 2007 
- January 2009 (Fall 2008). The upper limits on the 7-ray 
fluxes were considered as detections, with fluxes equal to 
one half of the limit. In "mid 2007" AGILE was pointed 
at 3C 454.3 when its optical main peak was already over; 
furthermore, we only have 5 7-ray points. The low statis- 
tics prevents us to obtain reliable results with the DCF 
for this period. In contrast, the period "fall 2007" offers 
a good opportunity to test the correlation, since the 7- 
ray flux, and even more the optical flux, exhibited strong 
variability. Moreover, the period of common monitoring 
lasted for more than a month. The corresponding DCF 
(Figure [T4|) shows a maximum DCF ~ 0.38 for a null 
time lag. 

However, the shape of the peak is asymmetric, and 
if we calculate the centroid ((Peterson et al.l IT998D . we 
find that the time lag is —0.42 days, i.e. about 10 hours. 
This result is in a g reeme nt with what was found by 
iDonnarumma et al.l (|2009l ) when analyzing the Decem- 
ber 2007 observations only, and implies that the 7-ray 
flux variations are delayed by few hours with respect to 
the optical ones. In the period "mid 2008" the main 
optical peak (and also the minor one) occurred when 
AGILE was not observing the source. Finally, we com- 



puted the DCF corresponding to the "fall 2008" pe- 
riod. We obtain a broad maximum, indicating a fair 
correlation (DCF max ~ 0.66), but with large errors, 
peaking at —2 day time lag, but with centroid around 
day. This result is consistent with that obtained in 
the "fall 2007" period, which appears to be the most 
robust one. Hence, for this case we estimated the uncer- 
tainty on the time lag by means of the statistical method 
known as "flux randomization / random subset selec tion" 
(FR/RSS [Peterson et all [1991 IRaiteri et al.1 12003ft . We 
run 2000 FR/RSS realizations and for each of them calcu- 
lated the centroid corresponding to the maximum. The 
resulting centroid distribution shown in Figure [T4l allows 
us to conclude that the 7-optical correlation occurs with 
a time lag of r — — O^lgJ, the uncertainty correspond- 
ing to a lcr error for a normal distribution. This result is 
consistent with a recent analysis o f the public Fermi dat a 
and the optical SMARTS data bv lBonning et al.l (|2009f l. 

8.2. Variability analysis 

The observed variance of a light curve for a specific 
detector can be written as 

1 N 

i=l 

where x is the average value of the Xi measurements. 
Moreover, since we deal with different detectors, in order 
to take into account the different count rates in different 
energy bands, and to compare their variance, we con- 
sider the normalized variance, S 2 /x 2 . In order to com- 
pute the intrinsic variance of a source light curve, the 
measurement errors must be taken into account, since 
they contribute an additional term to the varianc e . This 
approa ch was treated in det ail by iNandra et al.l (1997) 
and by lEdelson et al.l (|2002[ ). who introduced the term 
of "excess variance" : 

a 2 xs = S 2 -a 2 , (6) 
where a 2 is the mean squared error, 

i=l 

and Oi are the measurement uncertainties of a light curve 
points Xi. 

Thus, the normalized excess variance, 
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Fig. 12.— 



3C 454.3 light curves at different energies (see Section |7. II for details) over the whole time period. 
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Fig. 13. — Comparison between the light curves in different 
bands. Panel (a), (b), and (c) show the light curves in the op- 
tical, millimeter, and 7-ray energy bands, respectively. 
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Fig. 14. — Discrete correlation function between the 7-ray and 
optical fluxes during the "fall 2007" period. The uncertainty in the 
time-lag can be computed according to the FR/FSS method. The 
inset shows the resulting centroid distribution (see text for details). 



can be used to compare variances between different ob- 
servations. 

In order to quantify the flux variability in different en- 
ergy bands, we computed the fractional root mean square 
(rms) variability amplitude, F var , defined as 



l S 2 



(9) 



(see also lVaughan et al"]l2003L and references therein). 

Figure [15] shows the fractional rms variability ampli- 
tude at different frequencies. The optical R band is the 
one showing the highest degree of variability. This is 
partly due to the higher sampling of the optical data 
with respect to the the other frequencies. Nevertheless, 



Fig. 15. — Fractional rms variability amplitude as a function of 
frequency. 



a possible trend (higher fractional variability amplitude 
at higher frequency) is also present. This possible trend 
of the fractional rms variability amplitude with the loga- 
rithm of the frequency w as observed in othe r sources too 
(see e.g. PKS 2155-304. [Zhang et al.ll2005l ). and it was 
interpreted as signature of spectral variability. A more 
systematic study of the variability properties of 3C 454.3 
will be addressed in a forthcoming paper. 

8.3. Radio VLBI vs j-ray Data 

Figure [T^] clearly shows a strong enhancement of the 
radio core flux starting about on MJD 54500. The high- 
est flux density is on MJD 54742 at 15 GHz and on 
MJD 54719 at 43 GHz. This variability is not well cor- 
related with the variability at higher frequencies: optical 
and 7-ray data data show more different flares in the pe- 
riod MJD 54400-54800 (see Figure Q3] panels (a) and (c), 
respectively). Moreover, the radio flux density increase 
is smooth and longer in time, while 7-ray and optical 
flares are evolving faster. 

At 230 GHz the flux density variability mimics the 
VLBI radio core properties to MJD 54600, when a large 
flux density increase is visible, with a peak at about 
MJD 54630. At this frequency the source remains in 
an active phase up to MJD 54700 (Figure [T31 panel (b)). 

This poses an interesting question as to the nature of 
such an increase of th e core radio flux. As reported in 
iGhisellini et al.l ()2007f ) it is likely that the emitting region 
is more compact and has a smaller bulk Lorentz factor 
closer to the supermassive black hole. We can assume 
that in the region active at 43 GHz, in the quiescent 
state , the jet Lorentz factor is T ~ 10 ([Giovannini et al] 
|2001[ ). To obtain the flux density increase of the core at 
43 GHz (from ~ 5 Jy up to 25 Jy) the Doppler factor has 
to increase up to S ~ 30. Such an increase requires that 
the source is oriented at a small angle 9 with respect to 
the line of sight, since a large change in the jet velocity 
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will produce a small increase in the Doppler factor. A 
Doppler factor S = 30 can be obtained if 9 = 1.5° and 
T = 20, corresponding to a bulk velocity increase from 
0.9950 to 0.9987 (note that a larger orientation angle, 
e.g. 9 = 3° with the same increase in the jet velocity, will 
produce a small change in the Doppler factor <5, from 16 
to 19). 

The presence of one or more new jet components is 
not revealed in the high resolution VLBA images, even 
if the most recent VLBA images at 43 GHz suggest a jet 
expansion near to the radio core starting from MJD ~ 
54600. Because of different properties (multiple bursts at 
high frequency, a single peak in the radio band) it is not 
possible to correlate the radio peak with a single 7-ray 
or optical burst. We can speculate that a multiple source 
activity in the optical and 7-ray bands is integrated in 
the radio emitting region in a single event. This event 
(see Figure ITOf has a clear flux density peak on MJD 
~ 54720 and we can assume that 43 GHz is the self- 
absorption frequency at that epoch. This scenario is in 
agreem ent with the one discussed by iKrichbaum et al.l 
<|2008t ) (see the ir Figure 3 ). 

According to iMarscheil (|1983| ) , the self-absorption fre- 
quency, the source size, flux density, and the magnetic 
field are correlated as follows, 

B = 3.2 x 10" 5 x 9 4 v% S~ 2 5(1 + z)' 1 Gauss, (10) 

where B is the magnetic field in Gauss, 9 the angular 
size in mas (note that 9 = 1.8 X HPBW, where HPBW 
is the half power beam width), v m is the frequency (in 
GHz) of the maximum flux density S m (in Jy), and 6 is 
the Doppler factor, respectively. Moreover, we assume a 
particular value (a = 0.5) of spectral index in the opti- 
cally thin part of the synchrotron spectrum. 

Thus, we can use the radio VLBI data at 43 GHz to 
constrain the physical properties in the region where the 
source will start to be visible at at this frequency. The 
angular resolution in the jet direction of VLBA data at 
43 GHz is ~ .14 mas corresponding (as discussed in 
lMarscheH fl983) to 9 < 0.25 mas. Assuming 5 = 30, we 
obtain B < 0.5 Gauss. 

It is reasonable to assume that when the source is even 
smaller, the emission in the radio band is not visible be- 
ing self-absorbed, and that the local magnetic field is 
B < 0.5 Gauss when we start to detect the radio emis- 
sion. The size of this region should be smaller than 0.25 
mas (about 2 pc). 

8.4. Spectral analysis 

The correlation between the flux level and the spec- 
tral slope in the 7-ray energy band was extensively 
studied by means o f the analysis of the EGRET data. 
Nandikotk ur et all (|2007t l showed that the behavior of 
EGRET blazars is inhomogeneous. Figure [TBI shows the 
AGILE/GRID photon index as a function of the 7-ray 
flux at different epochs. A "harder-when-brighter" trend 
seems to be present in the long time scale AGILE data. 

Further long term observations of 3C 454.3 and of 
other bright 7-ray blazars at different flux levels with 
AGILE and Fermi will be crucial to assess this topic. 

Different emission mechanisms can be invoked to ex- 
plain the 7-ray emission. In the leptonic scenario, the 
low-frequency peak in the blazar SED is interpreted as 
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Fig. 16. — AGILE/GRID photon index as a function of the 7- 
ray flux above 100 MeV. Number beside each points represents the 
epochs listed in Table [T] 



synchrotron radiation from high-energy electrons in the 
relativistic jet, while the high-energy peak can be pro- 
duced by IC on different kinds of seed photons. In the 
synchrotron self-Compton [SSC ] model fGhiscllin i et al.1 
119851 iBloom fc Marscherl [l996h the seed photons come 
from the jet itself. Alternatively, the seed photons can 
be those of the accretion disk [exte rnal Compton s catter - 
ing of direct disk radiation, ECD, iDermer et al.l (1992)] 
or those of the broad-line region (BLR) c louds [exter- 
nal Co mpton scattering from clouds, ECC. ISikora et akl 
(1994)]. The target seed photons can also be those 
produced by the dust torus surrounding the nucleus 
[external Compton sc a tterin g from IR-emitting dust, 
ERC(IR). ISikora et al.l ([20?)1 ]. 

We fit the SEDs for the different observing periods 
by means of a one-zone leptonic model, considering the 
contributions from SSC and from external seed pho- 
tons originating both from the accretion disk and from 
the BLR (detailed d escription of this model is given in 
iVittorini et al.l [2~009| ) . Indeed, emission from both of 
them were detecte d during faint states of the source 
(Rai teri et al.ll2007t ). 

The emission along the jet is assumed to be produced 
in a spherical blob with comoving radius R by accelerated 
electrons characterized by a comoving broken power law 
energy density distribution of the form, 



(7/7b) Ql + (7/7b) ah ' 



(11) 



where 7 is the electron Lorentz factor assumed to vary 
between 10 < 7 < 10 4 , a\ and ah are the pre- and post- 
break electron distribution spectral indices, respectively, 
and 7b is the break energy Lorentz factor. We assume 
that the blob contains an homogeneous and random mag- 
netic field B and that it moves with a bulk Lorentz Factor 
T at an angle Qq with respect to the line of sight. The rel- 
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Fig. 17.— 3C 454.3 SED centered on MJD 54617-54618. Black 
triangles, red (blue) squares, red (blue) circles, green circles, and 
black stars represent radio, MJD 54617 (54618) Swift/ WOT, 
MJD 54617 (54618) Swift/XKT, RXTE/PCA, and AGILE/GRID 
data, respectively. UV and X-ray data are de-reddened and cor- 
rected for Galactic extintion. The thin solid, dotted, dashed, dot- 
dashed, and the triple-dot-dashed, represent the accretion disk 
blackbody, the synchrotron, the SSC, the external Compton on 
the disk, and the external Compton on the BLR radiation, respec- 
tively. The thick solid line represent the sum of all the individual 
components. 

ativistic Doppler factor is then 6 — [T (1 — (3 cosOo)] -1 , 
where /3 is the usual blob bulk speed in units of the speed 
of light. 

Our modeling of the 3C 454.3 high-energy emission 
is based on an inverse Compton model with two main 
sources of external target photons: (1) an accretion disk 
characterized by a blackbody spectrum peaking in the 
UV with a bolometric luminosity for an IC-scattering 
blob at a distance rd = 4.6 x 10 16 cm from the central 
part of the disk; (2) a Broad Line Region with a spectrum 
peaking in the V band, placed at a distance from the blob 
of tblr = 4x 10 18 cm, and assumed to reprocess 10% of 
the irradiating continuum (|Tavecchio fc G hiscllinil [20081 : 
IRaiteri et al1l2007l l2008bh . 

These two regions contribute to the ECD and the ECC, 
respectively, and it is interesting to test the relative im- 
portance of the two components that can be emitted by 
the relativistic jet of 3C 454.3 under different conditions. 
We summarize here the main results of our best model 
for the different time periods. 

Table [5] shows the best-fit parameters of our modeling 
of SEDs corresponding to the following periods (see Fig- 
ure [[J: (SED1), MJD 54617-54618, when 3C 454.3 en- 
tered a phase of high 7-ray activity; (SED2) MJD 54673- 
54693, when the 7-ray flux was almost constant; (SED3) 
MJD 54800-54845, when the source flux was at the min- 
imum level (about 70 x 10 -8 photons cm~ 2 s _1 ). In Fig- 
ures [T71 HH1 and 033 the thin solid, dotted, dashed, dot- 
dashed, and the triple-dot-dashed, represent the accre- 
tion disk blackbody, the synchrotron, the SSC, the ex- 
ternal Compton on the disk, and the external Compton 
on the BLR radiation, respectively, while the thick solid 
line represents the sum of all the individual components. 
The insert of Figure [T9] shows the portion of the SED 
dominated by the contribution of the disk blackbody ra- 
diation, which clearly emerges since the source is a rela- 
tive low state. 



Fig. 18.— 3C 454.3 SED during the period MJD 54673-54693. 
Black triangles, multicolor squares, circles, and black stars rep- 
resent radio, Swift/ WOT, Swift/XKT, and AGILE/GRID data, 
respectively. UV and X-ray data are de-reddened and corrected 
for Galactic extintion. The thin solid, dotted, dashed, dot-dashed, 
and the triple-dot-dashed, represent the accretion disk blackbody, 
the synchrotron, the SSC, the external Compton on the disk, and 
the external Compton on the BLR radiation, respectively. The 
thick solid line represent the sum of all the individual components. 
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Fig. 19.— 3C 454.3 SED during the period MJD 54800-54845. 
Black trinagles, multicolor squares, circles, and black stars repre- 
sent radio, Swift/WOT, Swift/XKT, and AGILE/GRID data, re- 
spectively. UV and X-ray data are de-reddened and corrected for 
Galactic extintion. The thin solid, dotted, dashed, dot-dashed, and 
the triple-dot-dashed, represent the accretion disk blackbody, the 
synchrotron, the SSC, the external Compton on the disk, and the 
external Compton on the BLR radiation, respectively. The thick 
solid line represent the sum of all the individual components. The 
insert shows the portion of the SED dominated by the contribution 
of the disk blackbody radiation. 



We find that the three SEDs can be reproduced well 
by very similar parameters, the main difference being the 
shape of the electron distribution and the break energy 
Lorentz factor. We note that the observed minimum vari- 
ability time scale, of the order of half a day, is consistent 
with the minimum variability time scale (~ 10 hours) 
allowed by the model fit. Finally, we computed for the 
three different SEDs the total power carried in the jet, 
Pjct, defined as 

Pjet = L B + L p + L e + L Iad ergs" 1 , (12) 
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TABLE 5 

Input parameters for the mode l of SED1, 
SED2, and SED3. See Section [8751 for details. 



Parameter 


SED1 


SED2 


SED3 


Units 


a\ 


2.3 


2.5 


2.0 




C*h 


4.0 


4.0 


4.2 




7min 


30 


30 


18 




7b 


300 


280 


180 




K 


80 


80 


100 


cm -3 


R 


21.5 


21.5 


21.5 


10 15 cm 


B 


2 


2 


2 


G 


5 


34 


31 


34 




L A 


5 


5 


5 


10 46 ergs" 1 


id 


0.015 


0.015 


0.015 


pc 


e 


1.15 


1.15 


1.15 


degrees 


r 


20 


20 


20 




-Pjet 


3.2 


3.7 


2.5 


10 46 ergs" 1 



where L B , L p , L c , and L rad art; Llio power carried by 
the magnetic field, the cold protons, the relativistic elec- 
trons, and the produced radiation, respectively. We ob- 
tain a value of Pjot of 3.2x 10 46 ergs^\ 3.7 x lO^ergs^ 1 , 
and 2.5 x 10 46 ergs^ 1 for SED1, SED2 and SED3, respec- 
tively. The total power of the jet is lower at the end of 
the AGILE observing period, following the general trend 
of the multiwavelength light curves. 

8.5. Jet geometry 

The light curves in Figure [13] show a different behavior 
starting from the end of 2007 among the different energy 
bands. 

As shown in lVillata et al.l (|2009T) . a possible interpreta- 
tion arises in the framework of a change in orientation of 
a curved jet, yielding different alignment configurations 
within the jet itself. 

During 2007, the more pronounced fluxes and variabil- 
ity of the optical and 7-ray bands seem to favor the inner 
portion of the jet as the more beamed one. On the other 
hand, the dimming trend in the optical and in the 7-ray 
bands, the higher mm flux emission and its enhanced 
variability during 2008, seem to indicate that the more 
extended region of the jet became more aligned with re- 
spect to the observer line of sight. 

9. CONCLUSIONS 

The AGILE high-energy long-term monitoring of the 
blazar 3C 454.3 allowed us to organize a few multiwave- 
length campaigns, as well as monitoring programs at 
lower frequencies, over a time period of about eighteen 
months. Thus, we were able to investigate the spectral 
energy distributions over several decades in energy, to 
study the interplay between the 7-ray and the optical 
fluxes, and the physical properties of the jet producing 
the non-thermal radiation. 

The global view we obtained after one and a half years 
of observations can be summarized as follows: 

1. The 7-ray emission for energy E > 100 Me V is 
clearly highly variable, on time scales of the order 
of one day or even shorter, with prominent flares 
reaching, on a day time scale, the order of mag- 
nitude of the Vela pulsar emission, the brightest, 
persistent 7-ray source above 100 MeV. 

61 http://www.to . astro . it/blazars/webt/ 



2. Starting from October 2008, 3C 454.3 entered a 
prolonged mid- to low-level 7-ray phase, lasting 
several months. 

3. Emission in the optical range appears to be weakly 
correlated with that at 7-ray energies above 100 
MeV, with a lag (if present) of the 7-ray flux with 
respect to the optical one less than one day. 

4. While at almost all frequencies the flux shows a 
diminishing trend with time, the 15 GHz radio 
core flux increases, although no new jet component 
seems to be detected. 

5. The average 7-ray spectrum during the different 
observing campaigns seems to show an harder- 
when-brighter trend. 

6. Our results support the idea the dominant emis- 
sion mechanism in 7-ray energy band is the inverse 
Compton scattering of external photons from the 
BLR clouds scattering off the relativistic electrons 
in the jet. 

7. The different behavior of the light curves at differ- 
ent wavelengths could be interpreted by a changing 
of the jet geometry between 2007 and 2008. 

The simultaneous presence of two 7-ray satellites, AG- 
ILE and Fermi, the extremely prompt response of wide- 
band satellites such as Swift, and the long-term monitor- 
ing provided from the radio to the optical by the GASP- 
WEBT Consortium will assure the chance to investigate 
and study the physical properties of several blazars both 
at high and low emission states. 
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